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K*- . Abstract 

(N 

^ ! Proceeding from WKB quantization conditions, we derive a semiclassical expression 

for the current of fermions scattered off a propagating bubble wall in the presence of 

Q\ • longitudinal gauge field. It agrees with the expression used by Nasser and Turok in 

i-C I semiclassical analysis of instability of electroweak bubble walls with respect to longi- 

' ! tudinal Z condensation. We discuss the resulting dispersion relation for longitudinal 

1^ ■ Z field and show that light species are important for the analysis of stability, because 



of their large contribution to plasma frequency. 
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1 Introduction 

In a recent preprint |1|, Nasser and Turok have found that interactions of bubble 
walls propagating during a first order electroweak phase transition with fermions 
can lead to a CP violating instability on the walls, a condensate of longitudinal Z 
bosons. This result, if confirmed by comprehensive study, is very important because 
it bears upon the main question in the theory of electroweak baryogenesis, namely, 
whether successful electroweak baryogenesis requires physics beyond the standard 
model or it does not. Because of deviations from equilibrium that occur during the 
wall propagation, the fermionic current, which enters the equation of motion for the 
Z field and contributes to its effective mass, cannot be obtained from thermodynamic 
relations and should be computed by direct averaging of a microscopic expression. 

As a first approach to the problem, Nasser and Turok considered the WKB ap- 
proximation for fermions and the thin-wall limit. (They have also checked results 
against a fully quantum- mechanical solution.) The expression for the current density 
they used in the WKB case would be obviously correct if the local value of canonical 
momentum of fermion remained unchanged during adiabatic switching on of the Z 
field. For a non-uniform Higgs field (bubble wall), this is not so, hence, a justifica- 
tion for that expression is required. In this note we shall show how that expression 
can be derived from WKB quantization conditions for fermions via certain manipu- 
lations with partial derivatives, similar to those used in thermodynamics. We shall 
also discuss the resulting dispersion relation for longitudinal Z field and show that, 
in realistic case, not only the top quark but also light species are important for the 
analysis of stability, because of their large contribution to plasma frequency. 

2 Derivation of the formula for the current 

The microscopic expression for the current density of a given species in the thin-wall 
case is 

J(x)=^/„j„(x), (1) 

n 

where j„ is the current density for the n-th single-particle mode and /„ is the cor- 
responding filling factor. The use of single-particle modes is specific of the thin-wall 
approximation that assumes that particles do not interact with each other while scat- 
tering off the wall. Though for a planar wall the non-trivial coordinate dependence is 



only that for the direction orthogonal to the wall, we keep the full three-dimensional 
coordinate in (|I]) because some of the expressions below are more general than that 
special situation. 

The partial current densities j„ are found by averaging the current density operator 
over eigenfunctions of the corresponding Dirac equation. If the Dirac hamiltonian is 
H , the current operator is 

^W = -SIM ■ '" 

where Z is some component of the field (in our case, that orthogonal to the wall). 
The eigenfunctions ipn satisfy 

Hll)n = Enlpn ■ (3) 

By a variant of the Feynman-Hellmann theorem (see for example ref. [0) 

J.(x)=Cj(x)V'n = -^. (4) 

Hence, the partial currents are related to adiabatic, fixed n variations of the energy 
eigenvalues with respect to the field. 

We can transform the fixed n variations to fixed energy variations via 

6En ( dEn\ ( 6n \ 



5Z(x) \ dn ) ^ \6Z{-x)J ^ 

which is derived by using jacobians (cf. ref.0) 

D{E,n) D{E,n) D{E,Z) 



(5) 



(6) 



D{Z,n) D{E,Z) D{Z,n) 

The resulting expression for the current is quite general (once the thin-wall limit 
is assumed) but it is especially useful in the WKB approximation, because in that 
case the quantum numbers have simple representation through WKB quantization 
conditions. 

Like ref. |ll], we assume the potential barrier created for particles by the wall and 
the Z field together to be monotonic. The WKB quantization conditions in the 
direction orthogonal to the wall (the only non-trivial ones) are 

rLz/2 

Pz{z)dz = 2-Knz , (7) 

-Lz/2 
r2*(n) 

2 / Pz{z)dz = 2Tmz + C (8) 

-L./2 



for transmitted and reflected fermions, respectively. Here Pz{z) is the longitudinal 
component of fermion momentum at a given spatial point, as determined from the 
classical dispersion law by conservation of single-particle energy and transverse mo- 
mentum, riz is a positive or negative integer. The classical turning point z* of reflected 
particles depends on state n = {nx,ny,nz). A constant (up to exponential in —L^ 
terms) phase shift C depends on the boundary condition for reflected particles at 
—Lz/2. The dispersion laws are listed in ref.[|l|. 

A single-particle mode can be specified by the value of momentum at one of the 
infinities (more precisely, at -^Lz/2), p^^. That incident momentum is related to the 
energy of the mode by the usual, field-independent expression. So, at given Px.oo, 
Py^ooi the variation at fixed energy is variation at fixed ^2,00- 

Next, we use ra^ instead of n in (^). Going over to one-dimensional variations we 
obtain 

5En 1 (dEn\ f Sn^ \ r] f OEA f dpz{z)\ 



5Z{z) LxLy\dnzj^\5Z{z)j^^^ {2Tc)LxLy ydn^ J ^ \dZ{z) J ^^^ 



(9) 



where r^ is 1 for transmitted states, for reflected states with z* < z, and 2 for 
reflected states with z* > z. Note that the variational derivative with respect to Z is 
converted into ordinary one. The derivatives are taken at fixed Px,oo, Py,oo, which is 
not indicated explicitly in (H). 

Using d^) and another transformation of derivatives 

(dpz{z)\ (dpz{z)\ fdpz^oo\ 



\ ^ ' / pz, 00 ^ ^ / Zj \ / pz 

we transform the current density as follows 

J{z) = I drixdnydnjnjniz) 

dpx,ocdpy,oo ^ , ^ fdEn\ (dpz{z)\ (dpz 



(10) 



(27r)3 ''''^'^'\dnJAdpz,oo A dZ 



Z \ •~'yz,oo / z \ " / p^ 



dpx,ocdpy,oo ^ , ; OE fdpz,oo\ (dpz{z)\ (dpz 



(27r)3 ~ dpz,oo \ driz ) z\ dp^^oo ) z\ ^^ 

dpJA\ (dE\ 
dpz,^)^[dz)^^ • 






(11) 



The last line is the expression used in ref.[|l|. We have thus shown how it can be 
derived from WKB quantization conditions. 



3 Dispersion relation and plasma frequency 

The formula for the current allows us to compute the response of the non-equilibrium 
plasma to longitudinal Z field. Because in the WKB calculation of the current both 
the wall profile and the gauge field were assumed slowly varying, we are actually 
considering the zero-momentum limit of the response. Accordingly, the relation (pJ]) 
between the current and the field is local in space. 

The linearized equation of motion — the dispersion relation for the Z field is 

diag(6c;2 - m\{z), iJ) - QjI{z) = . (12) 

This is a matrix equation because of the mixing between Z and photon field A due 
to plasma effects; u}p{z) is the non-equilibrium plasma frequency matrix. 



,2/.N / dJ{z)/dZ{z) dJ{z)/dA{z) 

'" ' ' dJ,rn{z)ldZ{z) dJem{z)/dA{z) 



Z=A=0 



where J now denotes the total, summed over all species, current coupled to Z and 
Jem is the total electromagnetic current. 

Non-equilibrium effects in the plasma frequency are important only for the heav- 
iest fermion, the top quark, which interacts effectively with the bubble wall. Evalua- 
tion of the integral ( pT]) shows [Q that for non-zero wall velocity, fermions produce a 
contribution to the current which has a square- root singularity in moo ~ m{z), where 
rrioo is the mass of a fermion in the broken phase at infinity and m{z) is its mass 
locally (so the singularity is behind the wall). This singular term corresponds to a 
negative contribution to the upper-left entry of the plasma frequency matrix, 

^l 

>ii ^"^ 477 cos2 Ow [m^ - m2 (z)] 1/2 

where aw = fi'^/(4vr) is the weak interaction constant, u is the wall velocity. One 
can see that the typical value of Pz{z) that give rise to the singular term ( |T4t) is of 
order [rrioo — m{z)Y^'^. For the WKB approximation to be applicable, this should be 
large compared to the inverse thickness of the wall given roughly by the Higgs mass 
rriH (all masses are those at the phase transition temperature). So, within the WKB 
domain, the singularity provides an enhancement factor of order 

moo/[m^ - m^{z)Y^^ < nioo/mH , (15) 



A^„n(^) - --r—-T7^T-3. „./„Mi/. ' (14) 



which can be considerable for top quarks. 

However, to decide if non-equihbrium contribution of top quarks to the dispersion 
relation, whether obtained by WKB or fully quantum-mechanical means, leads to 
an instability of the Z = state, we should consider it together with contributions 
of other particles.^] In particular, numerous light species (those with masses much 
smaller than temperature) give an essentially equilibrium contribution to Up. The 
equilibrium contribution of fermions and transverse W bosons to Up can be obtained 
simply by orthogonal transformation of the diagonal matrix of plasma frequencies of 
the SU{2) X U{1) basis. With A^^ light fermionic generations, light transverse W, and 
small wall velocity, this contribution is 



Is^ + |s^) + 2c^ Ng (l - |s2) sc + 2sc 
- Is') sc + 2sc^ [iNg + 2) s'c 

where c = cos 9w, s = sm9]y. For our estimates, we took A^„ = 3 but subtracted the 




top-quark contribution from (16), considered W bosons as light, and neglected the 
role of longitudinal W and Higgs bosons. 

In the WKB approximation, the quantity of interest is the magnitude of the non- 



equilibrium top-quark contribution to a>p, the most significant part of which is (|14|), 
corresponding to the onset of instability. It is determined approximately from 

where A, B, C are the entries of the equilibrium plasma frequency matrix of light 
species. (In a fully quantum-mechanical, non-WKB calculation, the relation between 
the current and the field becomes non-local |jl| and so does the criterion for instability.) 
Estimating the plasma frequency of light species as described above, we find from (|l^ 
that the onset of instability is at 

- AcjJ^i(^) = ml{z) + 0.25g'Ty cos^ Ow ■ (18) 

For realistic values of (poo/T, where 0oo is the expectation value of the Higgs field 
in the broken phase, the second term on the right-hand side is of the same order as 
m| in the broken phase, for example, for 0oo = 1-5T used in ref. IH, about 45% of it. 



^The author is grateful to M. Shaposhnikov for discussion that led to these considerations. 



We conclude that light species are in general important for the analysis of stability, 
because of their large contribution to the dispersion relation. 

The author is grateful to N. Turok for a lot of useful correspondence and sending 
parts of the longer paper prior to publication, and to T. Clark, G. Giuliani, S. Love 
and M. Shaposhnikov for discussions. This work was supported in part by the Alfred 
P. Sloan Foundation and in part by the U.S. Department of Energy. 
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